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Abstract 

Based on a single particle model which describes the time evolution of the 
wave function during tunneling across a one dimensional potential barrier 
we study the proton decay of ^*^®Pb from excited states with non-vanishing 
angular momentum 1. Several quantities of interest in this process like the 
decay rate A, the period of oscillation Toso th® transient time tij., the tunneling 
time ttun and the average value of the proton packet position are computed 

and compared with the WKB results. 
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I. INTRODUCTION 


The recent interest in nuclei beyond the proton drip line, prompted the study of proton 
emission phenomenon from the ground state of such nuclei [^. Like in the case of a decay, 
this process consists in the tunneling of the proton across a potential barrier. However for 
this kind of radioactivity it is vital to add the centrifugal potential to the Coulomb barrier 
since the majority of emitted protons are likely to decay from states with £ 7 ^ 0 . 

In last years a series of theoretical investigations, based on the quasi-classical approx¬ 
imation have been carried out in this field [^|-^. The proton half-lives of observed heavy 
proton emitters were calculated and compared with the experimental ones and in some cases 
a good agreement was found. 

In the Gamow approach the decay is treated as a stationary process, the penetrability 
being given by the ratio of probabilities of finding the quantum particle on each side of the 
potential barrier. In this image the dynamical aspects are neglected. However quantities 
like the tunneling time are important in decay processes, in fission or in fusion reactions. In 
order to include the time evolution of a wave packet propagating in a classically forbidden 
region one need to solve the time dependent Schrodinger equation (TDSE). There have been 
some attempts to incorporate the time-dependency in the study of the proton decay i by 
expanding the wave function into a term describing the bound-state of the prepared nucleus 
'4^o{A) and a set of exit channel wave functions x orthogonal to it which decomposed into the 
intrinsic states of the daughter nucleus — 1) and of the proton (l)p. However, imposing 
that the decay width is smaller than the resonance (quasi-stationary) energy (T -C Eq) the 
treatment reduce to the stationary Schrodinger equation for ‘ipp{r,t = 0) with a complex 
energy Eq — iT/2. Afterwards the usual computational procedure is carried on ^j. 

In this paper, based on a previous application of TDSE to the study of a decay [^H, 
we address the question of proton-decay from orbitals with £ 7 ^ 0 for the spherical nucleus 
208pb_ This process is different from the ground-state proton emission beyond the proton 
drip line. However the theoretical techniques that we employ in this paper can be applied 
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in both cases without major changes. In our time-dependent approach the description of 
the decaying system is fully contained in the state vector. The time evolution of this vector 
enables us to determine the decay probability at any moment t. We are interested to study 
the influence of the potential on the tunneling quantities and to determine to what extent 
the errors of the WKB method depends on the energy of the quasi-st at ionary state and on 
the angular momentum. 


II. DYNAMICAL APPROACH TO QUANTUM TUNNELING 


The interaction between the proton and the daughter nucleus jg (described by an 

average Woods-Saxon (WS) held which accounts for the nuclear potential 

w 


VN{r) = 


l + exp(^) 


( 1 ) 


a Coulomb potential, which is approximated by the interaction between the point proton 
and the uniformly charged spherical core of charge Z — 1 


Vc{r) = 






r 


,'r <R, 


c 


( 2 ) 


r > 


c 


and the centrifugal barrier 

The WS interaction is dehned by the nuclear radius Rq^ = with ro= 1.25 fm, the 

diffuseness a=0.7 fm, and the depth of the central potential Vo= 58 MeV 0. Note that in 
other papers, the depth of the potential is not hxed, its value being adjusted to reproduce 
the experimental energy of the quasi-bound state The Coulomb radius is given by 

R^ = ro(A - 1)1/3. 

The initial wave function of the proton was chosen to correspond to a quasi-stationary 
state of the potential V{r) = l//v(r) + Vc{r) + Vcf{r), with positive energy E > 0. The decay 
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width of such a metastable state can be calculated using an analytic method developed 
by Gurwitz et ah |]^. In this approach it is assumed that the proton occupies a bound 
eigenstate in the potential well V{r) + e{r) which represents a slight modihcation of V{r). 
Therefore we take 


i)p{r,t = 0) = 0, 


iV+e) 

nl 


(4) 


where, is an eigenstate of energy corresponding to the Hamiltonian 


Ho = -—A + V{r)+e{r) 

2/i 


(5) 


where fi is the reduced mass of the daughter-proton system. The modihcation e{r) reads 


e{r) = V{rraax) + {r - Tmax) tan^ - H(r), r > r^ax 
0 , T "E ^max 


( 6 ) 


where 6 gives the slope of the potential barrier beyond the point Tmax at which V (r) attain 
its maximum. 

In what follows we shall consider only the wave functions with the highest eigen¬ 
value bellow the barrier Vb = V{rmax)- 

In Table I we list the heights of the barriers (Vb), their locations (r^ax) and the selected 
eigenvalues {Eni) for a given value of the angular momentum 

The next step consists in the resolution of the time dependent Schrodinger equation: 


where 


d 


h 

H(r) = -A + V(r) 

2/r 


(7) 


( 8 ) 


A numerical procedure based on the iterated leap-frog method, provides the solution of this 
equation O- The equidistant spatial grid used in this method is large enough such that 
no interference of the outgoing and ingoing wave functions will take place during the time 
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intervals employed in the calculations. Once we get the wave-function which describes the 
time evolution of the proton packet through the potential barrier we are able to compute 
relevant quantities of the decay process. 

The tunneling probability can be expressed as the probability of hnding the proton 
beyond a certain point on the border which separates the zone inside the barrier from 
the external one 


POO 

PTD{rB,t)= |'0p(r,f)|V^dr 
Jtb 


(9) 


The decay rate is calculated according to the relation 

^TD{rB,t) = -—^— ( 10 ) 

It is also interesting to calculate the average value of the proton wave packet position 
operator, Vav, inside the nucleus 


rav{rB,t) 


r|^p(r,t)|V^dr 

\iljp{r,t)\^r^dr 


( 11 ) 


III. NUMERICAL RESULTS AND COMPARISON WITH WKB 

In Figure 1 a, b we present the time evolution of the proton wave function ‘ipp{r,t) 
for two angular momenta Z = 0, 5 at four different moments. We see the tendency of the 
wave function to decrease its amplitude in the interior of the barrier when time goes on. 
Eventually we observe that the fraction of the wave function which passed across the barrier 
behaves like a spreaded wave packet (at least on the spatial interval that we considered, i.e. 
up to 60 fm). This tendency is more pronounced for the wave function with i=5, because 
it faces a thiner barrier (see Table 2). It is worthwhile to mention that although the wave 
function amplitude decreases constantly in time its shape does not change to much inside 
the barrier. 

In hgures 2 a, b the time-dependent decay rate Xtd is ploted for angular momentum 
i = 0, 2, 5 and 8 . As has been pointed in an earlier work the decay rate undergoes 
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two regimes. In the first one, Xtd oscillates but increases on the average. The fact that 
at small times the decay rate is not constant as characteristic for exponential decay, but 
varies with time is typical for a quantum mechanical description [Q. This fact contrasts 
to the usual classical image which portrays the radioactive system as an ensamble of nuclei 
decaying independently one of each other with a probability which does not depend on time. 
In the second regime Xtd performs small fluctuations around a constant value, that we call 
asymptotic value Atd(oo). 

The decay rates presented in hgure 2 a correspond to tb = 11.6 fm, i.e. the inferior limit 
of integration is choosed to lay between the two turning points as can be observed from 
Table 2. If is increased, the irregularities occuring in the hrst regime are smoothed out. 
This fact is pictured in Figure 2 b, where tb is choosed to be 25 fm. 

In Figure 3 we represented the behaviour of the wave packet’s average position inside a 
potential region dehned by = 25 fm. As in the case of the decay rate we deal with two 
regimes. Whereas in the hrst regime Vav increases up to a certain limit, in the second one 
it performs small-amplitude huctuations around this limiting value as we expected since, 
according to a previous remark, the wave function does not change to much its shape. The 
period of oscillations in this second regime is denoted by Togc and its value can be deduced 
simply by measuring the distance between two maxima of Vav (right column in Fig.3). As 
has been noted in previous papers 0,§], the quasi-stationary state tends to penetrate the 
potential barrier by performing these small-amplitude oscillations instead of simply crossing 
the barrier from one side to the other. For this reason it makes sense to associate the 
frequency of collisions in the formula of Gamow with I/Tq^c instead of 1/Tcross, which is 
dehned by the quasi-classical expression (see eq.(13) bellow). 

In the WKB approximation the decay rate Xst is constant, its value resulting from the 
product of the barrier penetrability p and the collision frequency u 


XsT = iy-p 


( 12 ) 


where u is given by the inverse of the classical period of motion 
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( 13 ) 


2 /i dr 
h Jrti k{r) 


the wave number k{r) reading 


k{r) 


P(Q - r(r)) 


( 14 ) 


and the penetrability 

p = exp 1^-2^ dr^‘^{V{r) - Q)j (15) 

The stationary states computed in the modihed potential, become quasi-stationary when 
we turn on the real potential (without the modihcation e{r)) and their energy is no longer 
well dehned. Therefore, in all the above formulas the decay energy Q was taken to be 
the energy of the quasi-stationary state Eq, computed as the average energy of the time 
dependent Scrbdinger equation 


Eo = {'ipp{r,t)\H{r)\'ipp{r,t)) (16) 

In what concerns the decay rate, it can be inferred from Table 2 that the decay rates 
calculated in the time dependent approach X^d are in relatively good agreement with the 
WKB values \st{oo) especially when the quasi-stationary energy decreases with respect to 
the barrier height. Except the last case, with the highest angular momentum, in all other 
cases the WKB results overestimates the time-dependent values. 

The comparison between the times Tosc and Tcross (see Table 2) shows that they have 
very close values. However they are refering to different types of motion, the hrst describing 
the small oscillations of the proton wave function between the walls of the barrier during 
tunneling, the second, the classical movement inside the potential well. 

The transient time is dehned as the time interval between the moment when A^d starts 
to increase up to the moment when it reaches the limiting value Atd(oo). It depends on 
the energy of the quasi-stationary state. Its value can be deduced by inspecting hgures 2 
a and b. The tunneling time is related to the shift in time of the transition point between 
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the two above mentioned regimes, for the decay rate. Therefore it can be associated to the 
time necessary for the wave function to cross the barrier. Computing the decay rate for 
two different values of r^, i.e. the barrier’s turning points, and measuring the time delay 
between the two maxima of the two curves one gets the tunneling time (see Fig.4 a and b ). 
Notice that for £ = 0,2,5 the tunneling time decreases with increasing XtdIoo). However 
this does not happens for the state with higher angular momentum i = 8, where, although 
the decay rate is smaller, the tunneling seems to take place faster. In fact the tunneling 
time is correlated with the imaginary time, which is nothing else than the time necessary 
for the proton packet to cross the inverse potential, i.e. 



IV. SUMMARY 

Motivated by recent theoretical and experimental investigations on proton radioactivity 
phenomenon, we studied the time dependent characteristics of the proton tunneling from 
excited states in the spherical nucleus ^'^®Pb. Although this process is different from the 
ground-state proton emission beyond the proton drip line, the approach employed in this 
paper could be easily extended to the above mentioned case of proton radioactivity. 

Since other theoretical approaches are based on the semiclassical approximation we were 
interested to compare our exact results with the WKB ones. We found that the discrepancy 
between the two methods decreases when the difference between the top of the barrier and 
the energy of the quasi-stationary state increases. Our study does not concern a certain 
angular momentum state which could be measured in the decay reaction 208p|^*_j^207rjij _|_p 
Rather for a hxed set of Woods-Saxon parameters we investigated the dependence of the 
proton tunneling on time choosing one quasi-stationary state for every angular momentum. 
It seems that the accuracy of the WKB approximation increases with i. For a comparison 
with the experiment one should fit some of the WS parameters, e.g. the potential depth, in 
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such a way to reproduce the observed energy. However the present comparison between the 
WKB and the time-dependent approaches gives an idea of the error involved in the stationary 
approach and provides a good starting point for future investigations of proton decay using 
TDSE which could eventually answer to some questions related to this phenomenon. 
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TABLES 


TABLE I. The values of the barrier heights Vb-, their locations Vmax, the selected eigenvalues 
Eni and the wave functions number of nodes for different angular momenta L 


Angular Momentum 

Number of nodes 

rmax(fm 

) VB(MeV) E, 

.^(MeV) 

o 

II 


4 


10.62 

10.256 

7.78 

£ = 1 


3 


10.55 

10.631 

-0.04 

II 

to 


3 


10.45 

11.392 

7.23 

CO 

II 


2 


10.32 

12.560 

-2.06 

II 


2 


10.15 

14.163 

5.29 

£ = 5 


2 


10.00 

16.233 

12.53 

II 


1 


9.82 

18.800 

0.64 

II 


1 


9.62 

21.916 

8.44 

II 

00 


1 


9.45 

25.607 

16.58 

TABLE 11. The quasi-stationary energies Eq, 

the turning points of the potential (rti, rt 2 , rts), 

the asymptotic value of the decay rate Xtd{oo) and its WKB correspondent Xst, 

the oscillation 

period Tqsc and the crossing 

time 





Eo (MeV) ni (fm) 

rt 2 (fm) rts (fm) 

Xtd{oo 

) (s-') 

^ST (s 1) Tosc (s) 

Tcross (s) 

7.71 

9.24 

14.98 

1.25x10+20 

1.74x10+20 2.75x10-22 

2.95x10-22 

7.19 1.70 

8.90 

17.15 

1.90x10+10 

2.33x10+10 2.60x10-22 

2.74x10-22 

12.45 3.56 

8.74 

13.12 

2.23x10+20 

2.59x10+20 2.24x10-22 

2.43x10-22 

16.55 5.31 

7.89 

13.69 

2.79x10+10 

2.52x10+10 1.50x10-22 

1.79 X 10-22 
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Figure Captions 


Fig. 1. Time evolution of the squared wave function |'^p(r, for angular momentum 
(a) £ = 0 and (b) £ = 5 at moments t = 0, 1.5T0“^^, 3T0“^^ and 4.5T0“^^ seconds. In 
the second case the barrier through which the proton undergoes tunneling is thiner and 
therefore the probability to hnd it outside the barrier at a later time is larger. The wave 
function amplitude decreases in time and later on, the part of the wave function which 
already tunneled manifests itself as a well-spreaded wave packet. 

Fig. 2. The time dependent decay rates \td for the four selected quasi-stationary states 
of angular momentum ^ = 0,2, 5, 8. In the eq.(lO) we choosed (a) = 11.6 fm and (b) 

tb = 25 fm. In all four cases we observe that after a certain time the decay rate will fluctuate 
around an asymptotic value. In the second case these fluctuations are smoother. 

Fig. 3. The average value of the wave packet position operator Vav for three quasi- 
stationary states {i = 0, 2, 5). On the left side the asymptotic behaviour of Vav is observed, 
whereas on the right side we foccused on the small amplitude oscillations of the wave packet 
on its way to tunneling. 

Fig. 4. The decay rates for (a) i = 0,2 and (b) i = 5,8 when is choosed to be the 
internal turning point (full line) and the external turning point (dotted line). The difference 
between the maxima of the two curves gives the tunneling time ttun- 
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